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We study Andreev reflection in graphene nanoribbon/superconductor hybrid junctions. By using 
a tight-binding approach and the scattering formalism we show that finite-size effects lead to notable 
differences with respect to the bulk graphene case. At subgap voltages, conservation of pseudoparity, 
a quantum number characterizing the ribbon states, yields either a suppression of Andreev reflection 
when the ribbon has an even number of sites in the transverse direction or perfect Andreev reflection 
when the ribbon has an odd number of sites. In the former case the suppression of Andreev reflection 
induces an insulating behavior even when the junction is biased; electron conduction can however 
be restored by applying a gate voltage. 



PACS numbers: 74.45. +c, 73.23.-b, 72.10-d 

I. INTRODUCTION 

Graphene, a flat monolayer of Carbon atoms arranged 
in a two-dimensional (2D) honeycomb lattice, is a newly 
realized 2D electron system 1 which has attracted a great 
deal of interest because of the novel physics which it ex- 
hibits and because of its potential as a new material for 
electronic technology. One of the fascinating aspects of 
this system is that in most cases its potentialities for de- 
vice applications are intimately related with fundamen- 
tal aspects of quantum mechanics. Many of the trans- 
port properties of graphene which are at the heart of the 
design of new functional nanostructures originate from 
conservation laws of certain quantum numbers. For ex- 
ample, the fact that electrostatic barriers in graphene are 
perfectly transparent to electrons scattering with angles 
close to normal incidence (Klein-paradox) is explained 2 
in terms of conservation of pseudospin (the sublattice 
degree-of-freedom necessary to describe graphene's non- 
Bravais honeycomb lattice 1 ). 

Early investigations on transport properties of 
graphene have analyzed current and n oise in the pres- 
ence of normal leads. Recent studies^ 4 * 5 * 6 * 7 * 8 * 9 * have 
pointed out that novel interesting phenomena arise when 
graphene is interfaced to a superconductor (SC). In the 
seminal papers by Beenakker and coworkers^ it has been 
shown that the peculiar band structure of graphene 
gives rise to the appearance of specular Andreev reflec- 
tion (AR) 10 , a novel type of AR that is absent in ordinary 
metal/ SC interfaces. These studies paved the way to ex- 
perimental investigat ions of the proximity effectPE^ and 
of supercurrent flo w 13 ! 14 ! in graphene. On the theoretical 
side the results of Ref. 3 have been extended in a num- 
ber of ways, e.g. to graphene bilayers^ and to the case 
of interfaces with d-wave SCs^, whereas further studies 
have focused on the subgap structure of SC-graphene-SC 
junction^, on crossed AFP, and on magneto-transport 8 . 

Most of the theoretical analysis carried out so far de- 
scribe the graphene sheet as an infinite (or semi-infinite) 
2D plane, and identify two energy scales relevant for 
transport, namely the superconducting gap Aq and the 




FIG. 1: (Color online) A hybrid junction between a zigzag 
graphene nanoribbon and a superconductor. 



potential difference between the normal side and the su- 
perconducting side, which allows, for instance, to switch 
from the regime of ordinary AR to the one of specular 
AR. In graphene ribbons^, however, the finite size of the 
sample yields an additional energy scale characteriz- 
ing the mean energy spacing between the ribbon bands. 
Since the typical ribbon size varies from 10 nm up to 
1 /iiiP, S can range from 300 meV down to 3 meV, and 
it is thus larger than (or of the same order of) the typical 
superconducting gap Ao < 1 meV. As a consequence, S is 
expected to play an important role in electron transport 
and AR in graphene. 

In the present work we address this problem by ana- 
lyzing electronic transport through a hybrid junction be- 
tween a Graphene NanoRibbon (GNRJ^E^ and a gc, 
as sketched in Fig. [T] We consider the case of a GNR 
with zigzag edges, which has been shown to represent 
fairly well the behavior of an arbitrarily shaped edge^. 

Novel effects emerge due to the finite size of the GNR. 
We find that AR is strongly affected by the conservation 
of pseudoparity, a quantum number characterizing the 
GNR eigenstates, which depends on whether the num- 
ber Nw of sites along the transverse direction is even or 
odd. In particular, while in GNRs with odd Nw the AR 
coefficient is unity for subgap voltages, in GNRs with an 
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even Nw AR is totally suppressed. In the latter case AR 
can be restored to a finite value by applying a gate volt- 
age. Even-odd effects with the same physical origin have 
been also found in normal transport through graphene 
p-n junctions (valley- valve effect jM21\22\^ 

The contents of the paper are briefly described as fol- 



lows. In Sect.|II|we introduce the tight-binding Hamilto- 
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nian that we use to model the GNR and provide analyt- 



















ical expressions for the GNR eigenfunctions at arbitrary 
energy, briefly discussing how transport properties of the 
hybrid GNR/SC junction are calculated. In Sect. |III| we 
report and discuss our results. Finally, in Sect. |IV| we 
summarize our main conclusions. 



(a) N w even 



II. MODEL HAMILTONIAN, GNR 
EIGENFUNCTIONS, AND TRANSPORT 
PROPERTIES 

We model the hybrid GNR/SC system by means of 
a tight-binding approach. In order to interpret the re- 
sults presented below, we start by describing the eigen- 
states in the "normal" side of the junction (the GNR 
side). The GNR Hamiltonian TL is a matrix with ele- 
ments TLij = —jd(ij) + USij, where 7 c± 2.8 eV is the 
hopping energy between nearest-neighbor sites on 
the honeycomb lattice, and U accounts for a constant 
electrostatic energy, which can be controlled by a gate 
voltage. 

Denoting by a\ and CL2 the honeycomb-lattice basis, 
each unit cell is labeled by a vector r = n\a\ +7120*2, 
with rii and ri2 integers. The origin is located at the 
geometrical center of the GNR and Nw is the number 
of sites in the transverse direction e y . A suitable shape 
of the cells turns out to depend on whether Nw is even 
or odd, as depicted in Fig. [2] It can be shown that the 
(unnormalized) eigenfunctions of a zigzag GNR at a cell 
r = (x,y) are given by^l 

*a(b), s (s, y) = e ik * XA ^ x) <s> A (B)Ay) (!) 

where xa(b)( x ) is the coordinate of the A(B) site of the 
cell along the longitudinal direction e xi $>a(b),s is the 
wavefunction along the transverse direction at A(B), and 
s =|,| is a spin label. Explicitly, xa(b) = x if Nw is 
even, whereas xa(b) — x T aV3/4 if Nw is odd (a ~ 
1.42 A is the Carbon-Carbon distance). Furthermore 

(*A,s(v)\-{ sm[k y (W/2 + y)} \ 

\ *bAv) )~\V (W/2 - y)] ) W 

where W = 3a mt[(Nw + l)/2] approximately coincides 
with the width of the GNR, with mt[x] denoting the in- 
teger part of x. The eigenvalues related to ([!]) read 

e = U + r] 7 3?e f(k x aV3, 3k y a) (3) 

where 

/(£, k) = e iKNw/2 [l + 2e iK/2 cos (^/2)] , (4) 



e y (b) N w odd 




FIG. 2: Coordinates and elementary cells in GNRs with (a) 
an even number of sites in the transverse direction e y and (b) 
an odd number of sites in the e v direction. 



and k y is the solution of ^sm f(k x ay/3, 3k y a) = (see also 
Ref. [T8j) . The quantity k x (k y ) is the component of the 
momentum along e x (e y ). Importantly, beside the mo- 
mentum, the eigenstates defined in Eqs. ([l])-(|2| exhibit an 
additional quantum number 77 = zbl, which determines 
whether the transverse cell wavefunction (&a + 3>#)/2 
is even or odd with respect to the longitudinal axis of 
the GNR. One can therefore term 77 as pseudoparity. 
Note that, in contrast to what happens in bulk graphene, 
the states ([l])-([2| are not eigenstates of the pseudospin- 
projection operator a • h for any unit vector ft. 

The eigenvalues £ as functions of k x describe the band 
structure of the zigzag GNR, where two bands exhibit 
dispersionless zero modes and the other ones are some- 
what reminiscent of the bulk-graphene Dirac conefl, as 
illustrated in Fig.[3]for the case U = 0. Energies are mea- 
sured with respect to the Fermi level, and solid (dashed) 
lines refer to particles (holes) bands, which are degen- 
erate for this particular case. The value of the pseu- 
doparity 77, shown inside the circles, alternates from a 
particle band to the next one, and takes opposite values 
in particle and hole bands. Importantly, in the range 
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k x a^/?> G [0,7r], the pseudoparity of a GNR with even 
Nw is opposite to the one for the case with odd Nw 

The energy separation S between the first band and 
the zero-mode energy (Dirac level) is given by 5 ~ 
9ttj/(6N w - 4) for even N w > 1 or S ~ 3ttj/(2N w - 2) 
for odd Nw ^ 1 • A finite J7 shifts the energy of the 
Dirac level away from the Fermi level, as shown in Fig. |4j 
breaking particle-hole degeneracy. 

The presence of the superconducting electrode is ac- 
counted for by the Bogoliubov-de Gennes HamiltoniarPSl 



T^BdG 



H-e F A 
A* s F - TL* 



(5) 



where T~L is the particle Hamiltonian, and e F is the Fermi 
level of the GNR/SC junction at equilibrium. The SC (s- 
wave) order parameter, which couples particles and holes, 
is described by a matrix A with entries A^ = A^j. 
Here A^ is taken to vary smoothly across the junction, 
between its maximum value Ao in the bulk of the SC and 
zero in the bulk of the GNRP. 

Within the scattering approach^ the zero-temperature 
differential conductance is given by 28 



G(V) 



4e 2 



[N(V) - R(V) + R A (V)] 



(6) 



where V is the bias voltage applied across the junction. 
In Eq. ([6| the prefactor 4 is due to spin and valley degen- 
eracies, whereas N(V), R(V) and Ra(V) are the number 
of transverse propagating modes (open channels) avail- 
able at energy eV measured from the Fermi level £p, the 
normal reflection coefficient, and the AR coefficient, re- 
spectively. Below the gap (eV < Aq) only AR processes 
can contribute to the conductance, since the unitarity of 
the scattering matrix implies that N — R = Ra in this 
range. 



III. RESULTS AND DISCUSSION 

We present our numerical results focusing on the phys- 
ically realistic regime S ^> Ao- We first discuss the case 
where the Fermi level lies at the Dirac level (U = 0) and 
then the case with a finite gate voltage U ^ 0. 
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FIG. 3: (Color online) The band structure of a GNR as a 
function of the longitudinal wavevector k x . Only the three 
lowest bands are illustrated. Energies are in units of 7. The 
solid (dashed) lines represent the particle (hole) bands. In this 
figure [7 = (zero gate voltage). In this case particle and hole 
bands are degenerate (for the sake of clarity the hole bands 
have been slightly shifted to the right). The signs + and — 
inside the circles indicate the value 77 of pseudoparity for each 
band, (a) Nw is even. In the range < e < S no AR is 
possible, since left-moving hole states "2" have pseudoparity 
opposite to the one of incoming right-moving particle states 
"1". In contrast, for e > 5, left-moving hole states "4" and 
"5" with the same pseudoparity as incoming particle states 
"3" are available, and AR is finite, (b) Nw is odd. In this 
case AR is possible in the range < e < S since left-moving 
hole states "2" have the same pseudoparity of incoming right- 
moving particle states "1". For e > 5 left-moving hole states 
"4" and "6" with the same pseudoparity as incoming parti- 
cle states "3" are available and AR is again possible. Note 
that left-moving hole state "5" in this case has pseudoparity 
opposite to the one of incoming right-moving particle states 
"3". 



A. Zero gate voltage 

In Fig. [5|a) we plot G, Ra, and A/", as functions of the 
bias voltage applied across the junction, for a GNR with 
even Nw, and in Fig.[5|b) we plot the same quantities for 
a GNR with odd Nw In the subgap regime eV < Ao a 
striking dependence on the number of sites arises. While 
for even Nw the AR coefficient Ra vanishes, the oppo- 
site occurs for odd Nw, where Ra is unity. Because for 
eV < Ao electron transport is possible only by virtue of 
AR processes, the conductance G also vanishes in this 
voltage range for even Nw- In contrast, for odd Nw, G 



takes its maximum value in the same voltage range. In 
the regime Aq < eV < J, where quasi-particle transmis- 
sion becomes possible, Ra drops abruptly to zero also 
for odd Nw This is due to the fact that the transmis- 
sion processes, which are characterized by a small mo- 
mentum transfer (intra- valley), largely dominate over AR 
processes, which instead involve large momentum trans- 
fer (inter- valley). We thus find a finite conductance with 
a value equal to the number N of open channels (in units 
of 4e 2 /h). Finally, for eV > J, Ra is finite but quite 
small in both cases (for clarity Ra has been multiplied 
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FIG. 4: (Color online) Same as in Fig.|2]but for U / (finite 
gate voltage). In this case for e < U incoming particle states 
and outgoing hole states have the same pseudoparity, and AR 
is allowed. 
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by an enhancement factor), since voltages in this range 
are well above the superconducting gap. 

These features (Ra = for eV < 5 if Nw is even and 
Ra = 1 for eV < Ao if Nw is odd) can be understood in 
terms of a pseudoparity superselection rule for the scat- 
tering states of the GNR. Let us consider a right-moving 
incoming electron (labeled by "1" in Fig. |2| with an en- 
ergy lying between s-p and the bottom of the second band. 
The only hole state available for an AR process is the one 
labeled by "2" , which is characterized by a pseudoparity 
r] that is (i) opposite to the one of the incoming electron 
for even Nw and (ii) equal to the one of the incoming 
electron for odd Nw 

A superconducting parameter that is constant along 
the transverse direction cannot couple states with op- 
posite pseudoparity (i.e. the pseudoparity of an electron 
impinging onto the SC interface cannot be flipped). Thus 
for even Nw the AR process is forbidden and Ra van- 
ishes, whereas for odd Nw this process is allowed and Ra 
is finite. Notice that in the latter case normal reflection is 
forbidden by pseudoparity conservation, yielding Ra = 1 
in the subgap regime. For A < eV < S normal trans- 
mission is the dominating process, so that Ra is strongly 
suppressed for both even and odd Nw A finite AR is 
restored for eV > S. Indeed in this case intra-valley scat- 
tering into hole states with the appropriate pseudoparity 
["4" in Figs. |3ja) and Fig. ^b)} is available. This su- 
perselection rule is the ribbon counterpart of the bulk 
Klein-paradox selection rules based on pseudospin 2 ^ and 



FIG. 5: (Color online) Transport properties of a GNR/SC 
interface. The differential conductance G in units of 
4e 2 /h (solid line), the total Andreev-reflection coefficient Ra 
(dashed line), and the number N of open channels (dotted 
line) are shown as functions of the applied bias eV (in units 
of Ao). Note that the horizontal axis has been broken in two 
ranges. The GNR has a width of about 50 nm, which cor- 
responds to 5 rsj 52 meV = 52 Ao (the superconducting gap 
has been fixed at the value Ao = 1 meV). (a) Nw = 250. 
The AR coefficient is non-zero only for eV > 5, and has been 
multiplied by a factor 500 for clarity, (b) Nw = 251. The 
AR coefficient is unity for eV < Ao and drops abruptly to 
zero for eV > Ao, remaining zero up to eV = S. The AR 
coefficient is again finite for eV > 8, and has been multiplied 
by a factor 200 for clarity. 



also explains the valley- valve effect in p-n junctions in 



B. Finite gate voltage 

Let us now consider the case U ^ (and U <i). In 
Fig. [6] we plot our results for G, Ra, and TV. 

A comparison between Fig. [5]^a) and Fig. [6^a) shows 
that the application of a gate voltage has dramatic con- 
sequences on electron conduction in GNRs with even 
Nw- This is evident from Fig. |4^a). AR processes are 
allowed by pseudoparity conservation in the bias range 
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in Fig. [5jb) , due to a slight momentum mismatch in the 
transmission channel between GNR and S sides, arising 
when U ^ 0. 

Before concluding we would like to comment on the 
double-spike structure of Ra at voltages eV = S ± U, 
which is seen in both panels of Fig. [6| This is due to the 
opening of new channels for AR processes made avail- 
able by the existence of the first-excited energy band (see 
Fig.|§. 



IV. CONCLUSIONS 



2 



1 1 ' — 1 ' 1 1 1 — 








(b) N w odd 








G 








\ 










|; /? A x5oo 







J — //— 







eV/A Q 



FIG. 6: (Color online) Same as in Fig.[5]but for U = 0.5 A (in 
this figure we have not plotted the number of open channels) . 

(a) Nw even. In this case the AR coefficient is restored to a 
finite value for eV < U, while it remains zero for U < eV < 5. 

(b) Nw odd. In this case no big qualitative changes are seen 
to occur by switching on a finite value of U. 



eV < U, while they are forbidden for U < eV < 5. 
This fact is reflected in the results shown in Fig. [6ja), 
where the differential conductance is finite (and equal 
to 8e 2 /h) for eV < U and vanishes in the range^ 
U < eV < \]U 2 + Aq. For subgap voltages electron con- 
duction can be switched on and off by the application of 
a gate voltage. The junction between an even-Nw GNR 
and a SC exhibits the operational behavior of an electron 
transistor based on the pseudoparity conservation law. 

In contrast, for odd Nw the application of a gate volt- 
age is expected to have no major qualitative impact on 
the transport properties of a junction between a GNR 
with odd Nw and a SC. Once again, this expectation can 
be traced back to the pseudoparity quantum number, as 
depicted in Fig. |4jb). The data reported in Fig. [6^b) 
show indeed the same qualitative behavior reported in 
Fig. [5jb). The only slight modification is that the drop 
of Ra to zero, in the range eV > Aq, is less abrupt than 



In summary, we have studied Andreev reflection in 
graphene nanoribbon/superconductor hybrid junctions. 

We have reported analytical expressions for the eigen- 
functions of the tight-binding Hamiltonian describing the 
graphene nanoribbon, which are valid at arbitrary energy 
and carry explicitly a definite pseudoparity. The super- 
selection rule stemming from this quantum number has 
two main implications on transport through the hydrid 
junction. For narrow nanoribbons with an even num- 
ber of sites in the transverse direction we have found a 
complete suppression of Andreev reflection in a range of 
energies that is huge when the Fermi energy lies at the 
Dirac level. This implies zero conductance at subgap 
voltages, which can however be restored applying a finite 
gate potential, opening up potential technological appli- 
cations of these hybrid junctions as electron transistors 
as well as nanorefrigerators. In contrast, in the case of 
narrow nanoribbons with an odd number of sites we have 
found perfect Andreev reflection at subgap voltages and 
an abrupt suppression of it at supergap voltages. 

The study of non-ideal edges, realistic interfaces, and 
bulk disorder is beyond the scope of the present work 
and is postponed to a future publication^. The role 
of electron-electron interactions and/or next nearest- 
neighbor hopping, which also has not been addressed in 
the present work, can be qualitatively understood along 
the following lines. These effects have been show n to lead 
to the opening of a gap A g at the Dirac leveP^, which is 
however typically much smaller than 5. Our conclusions 
thus remain valid for a large range of energies even when 
these effects are taken into account. 
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